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Abstract 

Probably the most natural energy functional to be considered for 
knotted strings is the one given by the electrostatic repulsion. In 
the absence of counter-charges, a charged, knotted string evolving 
along the energy gradient of electrostatic repulsion would progressively 
tighten its knotted domain into a point on a perfectly circular string. 
However, in the presence of charge screening the self-repelling knotted 
strings can get stabilised. It is known that energy functionals in which 
repulsive forces between repelling charges grow inversely proportional 
to the third or higher power of their relative distance stabilise self- 
repelling knots. Especially interesting is the case of the third power 
since the repulsive energy becomes scale invariant and does not change 
upon Mobius transformations (reflections in spheres) of knotted trajec- 
tories. We observe here that knots minimising their repulsive Mobius 
energy show quantization of the energy and writhe (measure of chiral- 
ity) within several tested families of knots. 



1 Introduction 



Knot theory was given a strong impetus when in 1860's Kelvin proposed 
that knots made out of vortex hnes of ether constitute elementary parti- 
cles of matter which at this time were thought to be atoms |^. However, 
development of atomic physics, first with classical and then with quantum 
models of atoms, failed to show a connection between knots and atoms. 
More recently though, today's elementary particles are again considered to 
be string like objects that may be closed ||2|, ^ and perhaps knotted. This 
string theory approach partially revives the ideas of Kelvin and provides a 
motivation for exploring quantization of energy and of other values in such 
physical systems such as knotted magnetic flux lines [Q, ||] or knotted soli- 
tons 



Knots made of self-repelling, infinitely thin strings in 3-dimensional Eu- 
clidean space have been considered by many authors j^, |8|, |^, |l^, 11, 



Usually it was assumed that knots were perfectly flexible for bending but not 
extensible. The repelling charge was assumed to be continually spread along 
the knots so that there were no elementary point charges and that the charge 
contained within a given curvilinear segment is proportional to the length 
of this segment. When the Coulomb energy would govern the evolution of a 
knot of this type one should observe that the knotted domain gets progres- 
sively tightened into a singular point while the rest of the knot would form 
a perfect circle. Progressive tightening of knotted domains in electrostatic 
knots may seem counterintuitive when one considers that the electrostatic 
repulsion grows inversely proportional to the square of the distance between 
every pair of charges. However, at the same time as tightening progresses, 
the length of the knotted domain decreases and therefore there is less and 
less charge in the tightly knotted domain. In fact the decrease of charge in 
the knotted domain is more rapid than its contribution to overall repulsion. 
Therefore the tightening that is driven by the decrease of the repulsion out- 
side of the knotted domain can progress until a knotted domain shrinks to a 
singular point and the entire string is perfectly circular. For this reason the 



Coulomb energy is not interesting as energy functional for knots [13 



However, it had been demonstrated that if the repulsion force were growing 
inversely proportional to the third or higher power of the distance between 
the repelling elements then the knotted domains in prime knots would have 
no tendency to shrink to a singular point |8|, |lO|, 12 1. The third power case 
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is especially interesting from mathematical and physical points of view since 
the energy of a knot becomes conformal and therefore does not change when 
the trajectory of the knot is rescaled or undergoes a Mobius transformation 
(reflection in a sphere) 14, In this work we study configurations of 



knots that minimise their Mobius repulsive energy. From now on we will call 
these configurations as Mobius knots. Several earlier studies used numer- 
ical simulations or analytical approach to investigate various properties of 
Mobius knots |^, 15|, however the relations between such characteristic prop- 
erties of Mobius knots as their energy, crossing number, writhe or average 
crossing number were not systematically examined before. 



2 Energy of Mobius knots 



Let us consider an unknotted closed string that has a few repulsive point 
charges equidistantly separated. An energy minimising shape of such an 
unknotted string would be then an equilateral polygon with the number of 
vertices corresponding to the number of the point charges. To be able to 
operate with a model of self-repulsive knot whose shape is independent of 
the number of charges in the knot one needs to assume that charges are 
not localised but continuously spread over the knot. This mathematical 
operation ensures that unknotted energy-minimising strings would always 
form a perfect circle independently of the level of carried charge. However, 
this non-physical assumption of continuous charge redistribution causes the 
energy of a knot to become infinite due to the repulsion of nearby elements. 
In order to correct for this problem of infinite energies O'Hara ||8| introduced 
a regularisation term and defined the energy 

EiK) = ff (-^ - dsK{x)dsK{y) (1) 

JJKxK \\x — y\ aK\X,y) / 

where dxix^y) is the shorter arclength distance within K from x to y. 
Notice that the integral of the second term corresponds to the repulsive 
energy of a straight segment with the same length and carrying the same 
charge as K. Another, computationally more stable approach is to neglect 
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tangential contributions to repulsion as nearest neighbour regions in smooth 
trajectories are practically co-linear and define the cosine energy 



E{K)=ff ^-L^^dsKix)dsKiy) (2) 
where a is the conformal angle between the tangents at points x and y. In 



fact it was demonstrated by Doyle and Schramm [15| that 



E{K) = E{K) - 4. (3) 

We have applied here the second approach and used Kenneth Brakke's 
program Evolver to obtain Mobius energy minimising configurations of var- 



ious knots and to calculate their energies (Ref. 16 but see also |http:/ /www 



susqu.edu/facstaff/b/brakke/). Examples of configurations of various knots 
minimising their Mobius energy can be seen in Reference 15. It should be 
recalled that actual shapes of Mobius energy minimizers of a given knot can 
substantially vary because all configurations obtained by Mobius transfor- 
mations from one Mobius energy minimizer are also Mobius energy mini- 
mizers. In practice these shapes depend on an arbitrary choice of starting 
configurations used for the energy descent. However, the actual energy val- 
ues obtained in our simulations converge to the same values independently of 
the starting configurations of a given knot. In addition we have checked that 
the obtained by us values for ten different {2,p) torus knots were at most 
different by 0.1 % from the values calculated using an analytical approach 
that can be applied to this class of knots ^ . 

Figure 1 shows the relation between the Mobius energy and the topolog- 
ical, minimal crossing number for knots belonging to six different families 
of knots. We have analyzed torus knots with Alexander-Briggs notation 
3i,5i,7i etc, twist knots with even number of crossings (Alexander-Briggs 
notation 4i, 6i, 8i etc), twist knots with odd number of crossings (Alexander- 
Briggs notation 3i, 52, 72, 92 etc) and three Conway families of knots with 
Conway notations (2p + 1, 1, 2p), {2p — l, 2, 2p) and (1, 2p — l, 1, 2p — l) where 
p are consecutive natural numbers (see Fig. lb for schematic explanation 
how these families are formed). Standard representations of these knots 
classified according to Alexander-Briggs notations can be seen in tables of 



knots |17, 18| while the Mobius energy minimising configurations are shown 
in Ref. 15 and one representative example of energy minimising configura- 
tions of torus knot lli is shown here in Figure Ic. From the data points 
in Figure la we have excluded founders of the families as these frequently 
belong to different families at the same time. It is visible that in all these 
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Crossing nuiriljer 



Figure 1: Energy of Mobius knots, a, Different families of torus knots, twist 
knots and Conway knots show the same slope. These analysed families are 
represented in b with k = 1 and p = 2. c. The difference of energy between 
sequential torus knots tends asymptotically towards 52.8 which corresponds 
to 26.4 for each new crossing. The inset shows one of the configurations of 
the torus knot lli that minimises the Mobius repulsive energy, its position 
on the plot is indicated. 
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knot families the energy grows with a practically identical rate that seems 
to be linear. Linear fit over the tested range indicates that for each new 
crossing the energy grows by around 26-27 units. To give an estimation of 
energy units it is good to point out here that within the energy defined by 
the equation (|T|) the energetic costs of closing an open string into unknotted 



circle is exactly 4 units [14]. Notice that some of the plotted lines in Figure 
la practically coincide with each other while other seem to be vertically 
shifted by a constant value. 

A closer look at energy values for the {2,p) torus knot family (3i,5i,7i 
etc.), including the founder 3i knot, suggests that the energy difference be- 
tween consecutive torus knots shows an asymptotic convergence toward a 
constant value of circa 53. Since the energy minimisation by simulation 
with Evolver (or other programs) has its well known limitations in descend- 
ing toward a global minimum, especially for complex trajectories, we have 
used an analytical approach to generate torus knots configurations that are 



believed to be minimizers of the Mobius energy 15|. Figure Ic shows that 
the difference of energy between sequential torus knots tends asymptotically 
toward 52.8 which corresponds to 26.4 for each new crossing. 

On the basis of our simulation data (Fig. la) and analytical approach 
(Fig. Ic) we conjecture that within such families of knots that iteratively 
increase their interwound regions with double helix structure [15|, the dif- 



ferences in Mobius energy due to each new crossing tend to an universal 
constant value. 



3 Writhe of Mobius knots 



As already discussed, knots minimising Mobius energy do not have unique 
shapes since Mobius transformations can create infinitely many different 
configurations that minimise the Mobius energy for a given knot. How- 
ever, certain characteristic properties of curves in space are not changed 
by Mobius transformation. It was proven by Banchoff and White that the 
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absolute value of writhe for a given trajectory is invariant upon Mobius 



transformation |19]. Writhe (Wr) measures the extent of chirality of closed 
curves in space and therefore provides an interesting measure for knotted 
trajectories minimising a given energy functional. The writhe corresponds to 
the average difference between numbers of right- and left-handed crossings 
perceived when a given curve in space is observed from a random direction. 
The value of writhe (Wr) is usually calculated using the Gauss integral for- 
mula 

Wr = // ^J!^21lhi^dsK{x)dsKiy) (4) 

where u and v are the unit tangent vectors to X at a; and y, respectively. In 
the case of tight knots minimising their rope length and that are known as 



ideal, the writhe values showed a quantization [gO|, |2l|, and the quantum 
of writhe depended on the type of introduced crossings |23, |^]. We decided 
therefore to check whether knots minimising their Mobius energy also show 
a similar quantization of writhe. 

Figure 2a illustrates that Mobius torus and twist families of knots show a 
clear quantization of writhe where the writhe increase for torus knots is more 
rapid than for twist knots. While the differences of energy within a given 
family were showing an asymptotic descent toward a limiting constant value 
(Fig. Ic), the differences of writhe showed a specific constant value (within 
the accuracy of our computational approach) that seemed to be independent 
of the complexity of the knot. In case of ideal knots (see Fig. 2b) it was 
observed earlier that the slopes of writhe increase for torus knots and twist 



knots were described by simple relations: 1+x and 1—x respectively [23 



where x was the same for both families and seemed to correspond to 3/7 
[p3| , 24 1 and where the integer values 1 (or —1 depending on the handedness 
of crossing) is due to inter-coil crossing contribution while the noninteger 



value X is due to intra-coil contribution to writhe |23|. We observed here 
that the slopes of writhe increase for Mobius torus and twist knots also show 
opposing deflections from the slope of one, however the deflection value x is 
close to 0.2 instead of 0.43. Thus torus-type of turns with positive signs of 
crossings introduce a writhe of circa 1.2 per crossing. Twist-type of turns 
with positive signs of crossings introduce a writhe of circa 0.8 per crossing. 
Torus- and twist-type of turns can be recognised by the orientation along 
both arcs enclosing double-arc flelds in minimal crossing diagrams of a knot 
[ p3[ . Parallel and anti-parallel orientations characterise torus- and twist- 
type of turns, respectively [p3|. 



In the case of ideal knots it was observed that the writhe of achiral knots 
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Figure 2: Writhe of Mobius knots, a anrfb s/ioif i/ie quantization of absolute 
writhe of Mobius knots and ideal knots, respectively, the analysed families 
are torus knots and twist knots, c, Absolute writhe slopes for all tested by 
us families of knots. 
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Figure 3: Odd and even twist knots differ in the sign and type of crossing 
in the terminal clasp. Odd twist knots have terminal clasp with torus type 
of crossings orientations along the arcs enclosing this double-arc field is 
parallel) while even twist knots have twist type of crossings (orientations 
along the arcs enclosing this double-arc field is anti-parallel) Notice that 
sign of crossings in the terminal clasp also changes although signs and types 
of crossings in the inter-wound region remain the same. 



was essentially equal to zero [20, 24]. We observed here the same tendency 
for knots minimising their Mobius energy as all achiral knots tested by us 
(i.e. 4i, 63, 83, 89, 812, 817 and 813) had their writhe practically equal to zero. 



Figure 2c shows absolute writhe slopes for all families of knots tested 
by us. It is visible that within each of these families of knots there is ap- 
parently constant, specific increase of writhe as one analyses consecutive 
members of respective families. Interestingly the increase of writhe can be 
simply predicted by analysing the type of turns that are introduced while 
creating a next member of the family. So for example in the {2p + 1,1, 2p) 
Conway family, as one goes from a knot to its successor four new crossings 
are introduced: two are positive torus crossings which increase the writhe 
by about 2.4 and two are negative twist crossing which decrease the writhe 
by about 1.6. Therefore the predicted increase of writhe of about 0.8 when 
divided by four crossings gives us the observed slope about 0.2 (see Fig. 
2b). We similarly can predict and explain why knots from even and odd 
twist families of knots follow slopes with the same inclination but that are 
vertically shifted by exactly four units in relation to each other. Figure 
3 shows that change from an odd to an even twist family of knots impli- 
cates that two torus type of crossings in the terminal clasp are replaced 
by two twist type of crossings of opposite sign. Since the corresponding 
writhe contributions of these crossings are 1 + x and — 1 + x (or —1 — x 
and 1 — X for opposite handedness), where x is a constant, the resulting 
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absolute difference of writhe is two per crossing irrespectively of the actual 
value of intra-turn contribution (x) to writhe. Since there is a change of 
two crossings the global difference of writhe between twist even and twist 
odd families of knots results in the observed vertical shift of corresponding 
writhe slopes by exactly four. We can apply a very similar type of reasoning 
to explain why Conway families (2p + 1,1, 2p) and {2p — 1,1, 2p) differ in 
their corresponding values of writhe by exactly 2. Notice that the same type 
of explanation (where intra-turn contributions to writhe cancel) applies not 
only to Mobius knots but also to ideal knots. Figure 2b shows that in the 
case of ideal knots, odd and even twist knots also show the relative vertical 
shift by exactly four although the actual slope values are different. Axial 
trajectories of ideal knots (rope length minimizers) can be regarded as a 
limit of energy minimising configurations when the energy is taken to be an 
integral of an ever-increasing inverse-power of the radius of certain circles 
passing through three points of the curve |25|. Therefore this observation 
strongly suggest that minimizers of repulsive functionals with any exponent 
between 3 and oo will always have a constant shift of writhe slopes between 
certain families of knots. Therefore for example slopes of absolute writhe for 
even and odd twist knots will always show a relative vertical shift by exactly 
four (compare Fig. 2a and b) provided that compaired knots minimize the 
same repulsive functional with the exponent ranging between 3 and oo. 



4 Relations between energy and crossings 



There are two principal measures of crossings applied to knots. Minimal 
crossing number is a topological invariant and corresponds to the minimal 
number of crossings that any representation of this knot type can have in 
any orthogonal projection. Average crossing number (ACN) applies to a 
given rigid embedding of a knot and corresponds to the average number of 
perceived crossings (irrespective of their handedness) when this particular 
embedding is perceived from a random point on a sphere enclosing a given 
trajectory. The ACN value can be calculated using the unsigned Gauss 
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integral formula 

ACN=// (5) 

where u and v are the unit tangent vectors to -fC at x and y, respectively. 
Several studies considered theoretically the relation between the Mobius in- 
variant repulsive energy and the two measures of crossings mentioned above. 
Freedman, He and Wang demonstrated for example that the Mobius energy 
of a knot is at least 27r fold bigger than the minimal crossing number of this 



knot I IC ] . Figure 4a shows the relation between Mobius energy and minimal 
crossing number of all knots up to 9 crossings. It is visible that data points 
fill a "cone" and that for these relatively simple prime knots the lower lin- 
ear bound of the energy could be put at least at 7.477r times the minimal 
crossing number. One can also analyse the upper bound for the energy as 
the function of minimal crossing number. In figure la we have shown that 
within a given family the energy grows proportional to the crossing number, 
but as the crossing number increases, founders of new families arise and 
these can have very high energy as compared to the members of already 
established families. Thus for example a twist knot 8i has the energy of 
217.4 while the 8i8 knot founding a new family has the energy of 283.9. The 
upper bound of the energy does not follow a linear relation with the crossing 
number but approaches a 7/5 power law (see Fig. 4a). The ACN value is not 
invariant upon Mobius transformation, however, we have checked that ACN 
value is very robust and hardly changes upon multiple Mobius transforma- 
tions. Therefore it seems reasonable to investigate the relation between the 
energy and ACN of Mobius, knots. Freedman, He and Wang showed that 
the Mobius energy of a knot increased by 56/11 is at least 127r/ll fold times 
bigger than ACN of energy minimising configurations [|lO| . Existence of this 
linear lower bound demonstrates that the Mobius energy can not grow with 
a power lower than 1 as a function of ACN. It could though grow with the 
higher power. In fact our own data (Fig. 4b) indicate that data points for 
the relation between the energy and ACN fill much narrower "cone" than 
it was the case for the relation between the energy and minimal crossing 
number (see Fig. 4a). The lower bound of this cone can be described by a 
nearly linear function, this function is given by 

y = 22.25 • (6) 

whereas the upper bound is better described by a nearly 7/5 power law, the 
function is here 

y = 14.78 • x^-^^ (7) 
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Figure 4: Relation between the energy and crossings for all prime knots up 
to 9 crossings, a, The relation between Mobius energy and minimal crossing 
number for all prime knots up to 9 crossings. The upper on lower bounds 
were best fits of power law functions y = ax". Fits satisfied the condition that 
no experimental point was outside of the bounds, h, The relation between 
Mobius energy and ACN for all prime knots up to 9 crossings, c, The energy 
as a function of average crossing number within different families of knots. 
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This power law behaviour may seem inconsistent with a hnear growth of 
energy with the ACN within different analysed families of knots (Fig. 4c) . 
However as the ACN increases new founders of knots' families enter and 
they frequently start their families with a linear growth but at a higher 
level. Again 81 and Sig knots are good examples as the difference of their 
energy is much bigger than expected from a slightly higher ACN of Sig knot. 



5 Mobius knots and their relation to random knots 
and other knotted physical systems 



It was shown earlier that the average writhe for a population of random 
knots of a given type closely corresponds to the writhe of ideal knot of the 
corresponding type [pO| , p6|| . We have shown here that with the exception 
of achiral knots the writhe values of Mobius knots do not correspond to 
that of ideal knots. Although within given families there is a linear relation 
between writhe of ideal and Mobius knots these relations are not universal 
and different families will be related by different linear relations. Thus for 
example the writhe of torus knots grows more quickly for ideal knots than 
for Mobius knots but the opposite is true for twist knots. Therefore we can 
conclude here that the writhe of Mobius knots is not related to the writhe 
of random knots of corresponding type. 



Time averaged writhe value of randomly fluctuating knots of a given 
type seems to be independent of the length of random chains forming a 
given knot. This is not the case of ACN (average crossing number) as its 
value progressively increases with the length of a random chain. However, 
it was observed that for relatively simple knots the differences between time 
averaged ACNs of randomly fluctuating knots of a given type and of un- 
knots of the same chain size closely corresponds to the ACN of ideal knots 
of the corresponding types |14|. However this does not apply to ACN values 
of Mobius knots since these are significantly smaller than ACN values of 
corresponding ideal knots. We conclude therefore that Mobius knots in con- 
trast to ideal knots are not good predictors of certain physical properties of 
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random knots of a given type such as knotted polymer molecules undergo- 
ing a random thermal motion. However for other physical knotted systems, 
Mobius knots may better approximate their behaviour than ideal knots. If 
we imagine a charged knotted string of dimensions comparable to an ef- 
fective screening radius at given conditions then all pair-wise interactions 
within such a knot should be repulsive. However, due to screening interac- 
tions caused by counterions the repulsion would not follow the Coulomb law 
but would decrease more rapidly with the separating distance, approaching 
perhaps a cubic root dependence of the distance. Short flexible polymeric 
molecules like single-stranded DNA can make knots of dimensions compara- 
ble to effective screening radius at specific ionic conditions, such knots may 
approach then the Mobius behaviour. Finally one can entertain a thought 
about string-like charged elementary particles (electrons, for example) sur- 
rounded by short-lived mixture of other charged particles and antiparticles 
generated from quantum fluctuation of vacuum. Electrons may minimise 
then an energy that resembles the Mobius energy described here. If on the 
way to relaxation a complex self-repelling knot could undergo from time to 
time a strand passage and progressively simplify its type that would provide 
a physical system with natural quantization of the energy. 
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